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Aims. Using XMM Newton data, we investigate the nature of the X-ray emission in the radio relic 1253+275 in the Coma cluster We 
determine the conditions of the cluster gas to check current models of relic formation, and we set constraints on the intracluster magnetic field. 
Methods. Both imaging and spectral analysis are performed, and the X-ray emission is compared with the radio emission. 
Results. We found that the emission is of thermal origin and is connected to the sub-group around NGC 4839. The best-fit gas temperature in 
the region of the relic and in its vicinity is in the range 2.8 - 4.0 keV, comparable to the temperature of the NGC 4839 sub-group. We do not 
detect any high temperature gas, resulting from a possible shock in the region of the Coma relic. We therefore suggest that the main source of 
energy for particles radiating in the radio relic is likely to be turbulence. From the X-ray data, we can also set a flux upper limit of 3.2 x 10"'^ 
erg cm"^ s"', in the 0.3 - 10 keV energy range, to the non-thermal emission in the relic region. This leads to a magnetic field B > 1.05 fiG. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: Coma - Intergalactic medium - X-rays: galaxies: clusters 



1. Introduction 

Some clusters of galaxies show diffuse extended radio emission 
in their outskirts, indicating the presence of relativistic particles 
and magnetic fields in the intracluster medium (ICM). These 
radio sources are commonly called radio relics. Their steep ra- 
dio spectra indicate that the radiating particles have short life- 
times (~ 10** yr), linked to synchrotron and inverse Compton 
(IC) radiation losses, and thus are reaccelerated by some mech- 
anism that acts with an efficiency comparable to the energy loss 
processes. 

Current theoretical models predict that relativistic particles 
radiating in radio relics are powered by energy dissipated in 
shock waves produced during the cluster formation. This pic- 
ture is supported by recent numerical simulations on cluster 
mergers (Ryu et al. 2003), in which shock waves are present 
at the cluster periphery. It has been suggested that the elec- 
tron acceleration required to produce the reUc emission results 
from Fermi-I diffusive shock acceleration of thermal ICM elec- 
trons (EnBlin et al. 1998) or by the adiabatic energization of 
relativistic electrons confined in fossil radio plasma that was 
released by a formerly active radio galaxy (EnBlin & Gopal- 
Krishna 2001, EnBlin & Briiggen 2002). 

To test theoretical models. X-ray data are needed to deter- 
mine the physics of the ICM in the radio relic regions. In this 
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paper we analyze XMM-Newton data of the Coma radio reUc 
1253-1-275. This source is the prototype of cluster relic sources 
(Giovannini et al. 1991 and references therein). It is located 
at the Coma cluster periphery in the SW direction, beyond the 
galaxy group associated with NGC 4839, at ~ 75' from the 
cluster center. This location coiTesponds to a projected distance 
of ~ 2. 1 Mpc. The structure of the Coma radio relic is elongated 
with the major axis of ~ 840 kpc and is roughly tangential to the 
Coma cluster radius. The integrated spectrum is fairly straight 
between 150 MHz and 5 GHz, with a spectral index a - 1.18 
(Thierbach et al. 2003). 

We performed an imaging and spectral analysis of archival 
XMM-Newton data, to derive the properties of this relic and to 
shed more light on its physical origin. This is the first detailed 
study of the X-ray properties of a radio relic. 

We assume the concordance cosmology with Ho=70 km s"' 
Mpc Q,„ = 0.3, and Qa - 0.7. At the redshift of the Coma 
cluster, z = 0.023, the angular scale of 1' corresponds to a 
linear size of 28 kpc. 

2. XMM-Newton observations and data treatment 

The region of the Coma radio relic was observed with XMM- 
Newton for a total of 22 ksec on June 10 2003. We retrieved the 
data from the archive and pipeline processed them using SAS 
version 6.5. We followed the guidelines in the XMM-Newton 
ABC-guide (Snowden et al. 2004) for the data preparation. 
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We screened the XMM-Newton observation for time intervals 
of high background emission (flares), which can contaminate 
XMM-Newton observations. We examined lightcurves of the 
observations in two energy bands: first in the 10-12 keV (12-14 
keV) energy band for MOS (pn) data (see also Majerowicz et 
al. 2004) and subsequently in the 0.3-10 keV energy band (see 
Nevalainen, Markevitch & Lumb 2005 for details on flares). 
After the flare selection, exposure times are of 17(8) ksec for 
the MOS (pn) cameras. 

We used the weight method for the vignetting correction, 
attributing a weight factor to each event (Arnaud et al. 2001). 
Since there is extended low surface brightness emission in the 
region, we need to correct for the background. For this we used 
blank sky observations compiled by A. Read (Read & Ponman 
2003), which we screened for flares in the same way as in the 
Coma relic observations ' . 

3. Results 

3.1. Imaging analysis 

We chose for the imaging analysis the energy band 0.3-2 keV, 
since it allows the best signal-to-noise ratio (see also Scharf 
2002). Figure^shows the filtered count-rate image of the MOS 
cameras in this energy band. Two XMM-Newton observations 
were merged together for this image: the data taken in the PV- 
phase of the sub-group around NGC 4839 (Neumann et al. 
2001) and the observations described here. The PV-phase data 
(latest pipeline processing) were treated in the same way as 
the Coma relic observations. Figure is vignetting-corrected 
and background-subtracted using the above-mentioned blank 
sky observations. As can be seen from the image, the X-ray 
emission around the sub-group extends into the radio relic re- 
gion and shows a strikingly similar outer boundary. It is clear 
that the X-ray emission in the relic region comes from the sub- 
group. It cannot originate from the Coma cluster itself, as seen 
when extrapolating the beta-model derived by Briel, Henry & 
Bohringer (1992). In fact, the Coma cluster emission at these 
radii is more then 20 times fainter than what is observed. 

3.2. Spectral analysis 

To check for the presence of a shock wave and of non-thermal 
emission, we selected two regions for the spectral analysis (Fig. 
|5J: one ellipse of 600 x 330 kpc^ in the relic region (1) and 
another ellipse of 3 10 x 190 kpc^ to the NE of the reUc towards 
the Coma cluster (2). The total number of source photons in 
flie 0.3-10 keV band in regions (1) and (2) is 13520 and 10620, 
respectively. In these regions, the source count-rate in the 0.3- 
10 keV energy band is similar to the background count-rate. 

Table \l\ shows the results of the spectral fitting analysis 
with XSPEC. For this study, we used the latest response matrix 
files avalable at ESA's XMM-Newton homepage and the on- 
axis auxiliary reponse files describing the effective area of each 

' We also tried the background files from Nevalainen, Markevitch & 
Lumb (2005); however, using these files after additional flare screen- 
ing we obtained fit results with values higher than those obtained 
using Read's background files. 
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Fig. 1. Color-scale X-ray emission of the NGC 4839 group in 
the 0.3-2 keV energy range obtained from the MOS cam- 
eras. The image, obtained from a mosaic of two pointings, has 
been smoothed with a gaussian filter with cr - 10", and further 
with a median filter. Dark blue areas are regions with no X-ray 
data. The green contours represent the X-ray brightness at 8, 
1.6, and 2.4 x 10"^ cts pixel"' (pixel=4.4"x 4.4"), which 
correspond to 2, 3, and 5 cr level over the background. The 
radio emission from the relic at 90 cm is shown by the black 
contours, whose levels are 2.5, 4, 8, and 16 mJy/beam (beam - 
55"xl25", RAxDEC). 




Fig. 2. Zoom of Fig ^ on the region around the radio reUc 
1253-1-275. The radio contours are identical to Fig^ The el- 
lipses indicate the regions used for the spectral analysis. 

camera as a function of photon energy. We chose the corre- 
sponding background regions in the same detector coordinates. 
There is an intensity variability in the instrumental background 
of XMM-Netwon. To coiTect for this effect, we compared the 
source and background count-rates in the 10-12 keV (12-14) 
keV energy band for the MOS (pn) cameras over each entire 
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Table 1. Spectral fit results in the relic region (1) and in the 

region NE of the relic (2). For the thermal component we use 
XSPEC's MEKAL-model with free metallicity and hydrogen 
column density (nH). 



Reg. 


model 


kT/a 
keV/ 


nH 
lO^Ocm-^ 


red.^ 


1 
1 

2 


th. 
th./ n-th. 
th. 


•'••'-0.4 

3.2!;;-5/2.2 

Q Q+0.7 

■'•■'-0.5 


<0.65 
<0.01 


1.16 
1.17 
1.41 



Note. Confidence intervals of the parameters are at 90% level. 



camera (see Majerowicz et al. 2004 for details). The ratio of 
the count-rates is used as the normalization factor of the back- 
ground to the source spectrum. The normahzation value in our 
case is about 1.2, which indicates that the instrumental back- 
ground in our observation is higher than the instrumental back- 
ground in the blank sky observations. The chosen background 
normahzation of 1 .2 provides the lowest results. Varying the 
normalization factor, the changes of the spectral fitting results 
are always within the error bars. 

We fit both a thermal model and a combined thermal/non- 
thermal model with a fixed power law index of 2.2, derived 
from the radio spectrum (see Sect. 4.2), to the relic region. 
Both fits give consistent results and point to a thermal origin of 
the X-ray emission, since the combined fit gives a non-thermal 
component contribution consistent with zero. A spectral fit with 
only a non-thermal component of or = 2.2 does not provide 
an acceptable fit, giving a higher reduced and a nH value 
of 9 X 10^° cm"^, i.e. an order of magnitude higher than the 
observed value of 9 x lO'^ cm"^ (Dickey & Lockman 1990). 
Moreover, a thermal-model fit applied to the region NE of the 
relic provides a gas temperature fully consistent with that of the 
relic region. The relatively high reduced value in this region 
can be explained by residuals in the instrumental background 
subtraction, which are amplified by the vignetting correction. 

The non-thermal emission in the relic region, if at all de- 
tectable with XMM-Newton, would be much lower than the 
thermal emission. We calculate the upper limit of this non- 
thermal component in two ways: a) spectral fitting, which gives 
a maximum flux in the 0.3 - 10 keV band of 2.5 x lO"'-' erg 
cm"^ s \ and b) imaging analysis. For the latter, since no ex- 
cess X-ray emission directly related to the relic is detected, 
we can calculate how much of this emission can be "hidden". 
Assuming Poisson statistics, we estimate that this is 3 .2 x 10 
erg cm"-^ s~\ which is in excellent agreement with the value 
obtained from the spectral analysis. We adopt as an upper limit 
of the non-thermal component the more conservative value, i.e. 
3.2 x 10-12 ergcm-2 s'^ 

4. Discussion 

The results of this paper can be summarized as follows: 
a) There is extended X-ray emission in and around the Coma 
relic region, which is coimected to the sub-group around 
NGC 4839. The group is on its first infall onto the cluster, 
with velocity ~ 1700 km s as suggested from optical data 
by CoUess & Duim (1996), and as further demonstrated by the 



effects of ram pressure stripping detected from XMM-Newton 
data that indicate the infall direction (Neumann et al. 2001). 

b) The best-fit gas temperatures in the region of the relic and in 
its vicinity (NE region) are very similar, around 3.3 keV. This 
value is lower than the 4.4 + 0.4 keV found in the main body 
of file NGC 4839 sub-group (Neumann et al. 2001), although 
the diff"erence is marginally significant. However, what is most 
remarkable is that we do not detect any high temperature gas, 
which could indicate the presence of a shock at the location of 
the Coma relic. 

c) The emission of the relic is largely dominated by thermal 
emission. This is deduced by the spectral fit analysis and is also 
supported by the following arguments: i) the X-ray emission is 
connected to the hot gas of the NGC 4839 sub-group, and ii) 
there are no distinct X-ray features connected to the radio relic 
itself. We can only set a flux upper limit of 3.2 x 10 erg cm"^ 

in the 0.3-10 keV band for the non-thermal component. 
A brief discussion of the impUcations of our results follows. 

4.1. Merging process and relic formation 

According to current models, the relativistic electrons radiat- 
ing in radio reUcs are expected to be accelerated by shocks 
originating from cluster mergers. In a fully ionized plasma, a 
shock with compression factor C and Mach number M can ac- 
celerate particles to a power law distribution in momentum, 
with slope 5 = (C -H 2)/(C - 1) = 2(M^ + l)/iM^ - 1) 
(Drury 1983). The spectral index a of the radio emission is 
thus related to the Mach number of the shock by the rela- 
tion: a - (M^ + 3)/2(M^ - 1). In this framework, the Coma 
rehc with spectral index a = 1.18 should originate from par- 
ticles accelerated by a shock of Mach number M ~ 2. In the 
case of a monoatomic gas, the ratio between the post-shock 
and pre-shock temperatures T2 and Ti is given by T2/T1 = 
(5M^ - 1)(M2 + 3)/i6M- (see, for example, Sarazin 2002). 
This implies that the temperature jump across a shock of Mach 
number M ~ 2 is 2.1, which is not observed in our data. 

The above calculations refer to the simplest case. More rig- 
orously, one should consider a) the eff'ect of electron ageing on 
the observed steep radio spectrum, b) the influence of magnetic 
field on the shock parameters, c) the case of re-acceleration of 
relativistic particles, leading to possibly higher Mach numbers. 

No temperature gradient is detected between the region of 
the rehc (1) and the pre-shock region (2). This cannot be due 
to projection effects. In fact the high resolution radio images 
of the relic show a sharper brightness decrease on the west- 
ern (outer) edge (Giovannini et al. 1991), and the infall path 
of the NGC 4839 group to the cluster in a two-body model 
lies at about 74° to the line of sight (Colless & Dunn 1996). 
Thus, a strong inclination of the shock with respect to the sky 
plane, which would spread the observed temperature jump over 
a larger region, should be ruled out. The identification of re- 
gion (2) as the pre-shock region relies on the assumption that 
the group is on its first passage into the cluster core. This 
is a widely accepted scenario (see beginning of Discussion). 
Nevertheless, the low temperature at the relic location is a sig- 
nificant result in itself, in particular because it is lower than 
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that of the NGC 4839 group, thus indicating the lack of high 
temperature shocked gas in the relic region. 

This result favors the hypothesis that turbulence may be 
the major mechanism responsible for the physical origin of the 
radio relic emission. During its infall onto the Coma cluster, 
the sub-group around NGC 4839 encounters a region of rel- 
ativistic particles connected to a magnetic field. The interac- 
tion between the ionized moving plasma and the magnetic field 
would imply energy transfer from the ICM to the relativistic 
particles. As for the origin of relativistic particles, EnBlin et 
al. (1998) suggested that they could be produced by the tailed 
radio galaxy NGC 4789, located to the SW of the relic. 

The energy for particle acceleration comes from the kinetic 
energy of the X-ray gas, which is subsequently slowed down. 
We see indications that the X-ray gas slows down: for example, 
the gas of the sub-group is lagging behind the galaxies of the 
sub-group (Neumann et al. 2001). This deceleration is partly 
due to ram pressure stripping, which acts more on the extended 
X-ray gas than on the galaxies. The total minimum energy in 
the radio relic, computed under equipartition (Thierbach et al. 
2003), is ~ 3.5 X 10^*^ erg, assuming that the radio relic volume 
can be approximated by an ellipsoid of 30'xlO'xlO', corre- 
sponding to ~ 3.45 X 10^ kpc^. The kinetic energy of the ICM 
of the sub-group, assuming a gas mass of 10'-^ M0 and a ve- 
locity difference of 1700 km sec"', is ~ 3 x 10^^ erg. Thus, the 
energy from the group motion can be transferred to the rela- 
tivistic electron population with a reasonably low efficiency. 

The head of the NGC 4839 group crossed the relic region 
~ 5 X 10^ yr ago. However, the relativistic electrons radiating at 
frequencies between 150 MHz and 1.5 GHz are only 1 2.5 x 
10** yr old, thus supporting the view that particle reacceleration 
has occurred more recently. 

4.2. Limit to the magnetic field 

X-ray emission of non-thermal origin is predicted to origi- 
nate in a synchrotron source, due to the scattering between 
the radio-emitting electrons and the microwave-background 
photons. Since synchrotron and IC emissions are produced 
by the same electron population, they share the same emis- 
sion spectral index a. This spectral index relates to the in- 
dex 6 of the power-law electron energy density distribution 
as 6 = 2a + I, and to the photon index of the IC emis- 
sion as Fx = a -H 1. Combining the standard formulae of the 
synchrotron and IC emission mechanisms (see, for example, 
the derivation in Govoni & Feretti 2004), a lower limit to the 
volume-averaged magnetic field can be obtained from the X- 
ray flux upper limit for the non-thermal component. 

Using a radio flux at 6 1 MHz of 6 1 1 mJy (Giovannini et al . 
1991) and the spectral index of 1.18, we obtain a magnetic field 
B > 1.05 ijG. This value can be compared with the estimates 
obtained by other arguments. First, the equipartition magnetic 
field derived by Thierbach et al. (2003) by integrating the en- 
ergy spectrum of particles over all energies above 300 MeV, 
and by assuming a filling factor = 1 and a proton-to-electron 
energy ratio k = 1 , is 0.67 fiG, scaled to our cosmology. Second, 
the expected value of the magnetic field at the distance of the 



relic, derived from the profile of the Coma cluster magnetic 
field presented by Brunetti et al. (2001), is of about 0.1 fiG. 
The profile is obtained in the framework of a model, which re- 
produces the observational properties of the Coma radio halo 
invoking electron reacceleration. Finally, the average magnetic 
field in the Coma cluster, obtained from the non-thermal hard 
X-ray emission detected with BeppoSAX is ~ 0.2 fiG (Fusco- 
Femiano et al. 2004). 

The lower limit derived here is almost two times larger 
than the equipartition value and is about an order of magni- 
tude larger than the value derived from the profile presented by 
Brunetti et al. (2001). This indicates that the magnetic field in 
this region of the Coma cluster has undergone some amplifi- 
cation. This could be related to the sub-group gas overdensity 
(Dolag et al. 2001). Another possibility is that the magnetic 
field could have been deposited in the ICM by formerly active 
radio galaxies, also responsible for the origin of relativistic par- 
ticles. Since this is a very small region compared to the whole 
cluster, a high magnetic field here is not in contrast with the 
low average value of the Coma cluster, which should still be 
dominated by the magnetic field decline with radius (see also 
the discussion in Fusco-Femiano 2004). 
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